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ABSTRACT 

The masses of stars including stellar remnants are almost exclusively known from bi- 
nary systems. In this work, we study gravitational microlensing of faint background 
galaxies by isolated neutron stars (pulsars) . We show that the resulting surface bright- 
ness distortions can be used to determine the masses of neutron star. Due to different 
evolutionary histories, isolated neutron stars may have different masses from those in 
binary systems, and thus provide unique insight into their equation of states under 
extreme physical conditions. We search for existing pulsar catalogues and find one 
promising pair of a nearby pulsar and a background galaxy. This method will become 
more practical for the next generation optical and radio surveys and telescopes. 

Key words: Gravitational lensing: micro - stars: neutron - stars: pulsars 
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1 INTRODUCTION 

From stellar evolution, nume rous (~ 10 s ) neutron star s are 
expected in the Galaxy (e.g. lArzoumanian et alj|2002h . So 
far only a tiny fraction of these have been discovered, mostly 
as pulsars. Even a small fraction have measured masses, usu- 
ally for binary systems (due to additional orbital informa- 
tion). However, the masses of these neutron stars may well 
be affected by physical processes (e.g. accretion) between the 
companions. The masses of isolated stellar remnants (such as 
neutron stars) are thus of particular interests since they may 
be more easily compared to that predicted by stellar evolu- 
tion; such comparison can provide important constraints on 
the equation of state under extreme physical (e.g. density) 
conditions. 

The recent discovery of a 2Mq neutron star in a 
binary system alread y rules out some equation of state 
(|Demorest et alj|2010h . More mass determinations of neu- 
tron stars in isolated or binary systems would be important 
for such constraints. Gravitational microlensing has already 
demonstrated its capabilitie s to measure the masses of iso 
lated normal stars 



e.g. 



Smith et al.ll2003; Gould et al 



120091 : iHwang et all bold : iBatista et alj|2009l ; see iMad 12012 

for a review). 

IPaczvnskil |l995l ) studied how the masses of nearby 
dwarfs with high p roper motio n s can be determined using 
microlensing, while rColes et all (|2010l ) studied how "weak" 
microlensing of background galaxies can be used to deter- 
mine the masses of nearby stars. For the usual "strong" mi- 
crolensing, we observe the magnification change as a func- 
tion of time; in "weak" microlensing, we observe the sur- 
face brightness distortions of the background galaxy by the 



nearb y lens as it moves across the galaxy. The lColes et al.l 
(|20ld ) study has one potential difficulty, that is, their nearby 
stars will be very bright compared to the background galaxy, 
and so some type of coronagraph has to be used to block the 
light from the foreground star. In this work, we study gravi- 
tational microlensing of background galaxies by nearby stel- 
lar remnants (< lkpc). The stellar remnants are expected 
to be faint in the optical, and so the "glaring" problem will 
be much less severe. 

The outline of this p aper is as follows. In [J2] we first 
briefly recall the method of lColes et all (|2010l ). In Sg] we per- 
form simulations for the Thirty Meter Telescope (TMlQ). In 
!j4] we perform a preliminary search for pulsars overlap with 
background galaxies, and present some preliminary matches. 
In 35] we discuss our results further. 



2 METHOD 

In principle if we have an image of the background galaxy 
at the time before the lens star crosses in front of the galaxy 
and a second image when the lens star is passing the galaxy, 
a simple subtraction of the two images will reveal the dis- 
tortions by the lens. From these distortions, we can derive 
the angular position of the lens star and the corresponding 
Einstein radius. If the lens distance is known, then the mass 
of the lens star can be directly estimated from the derived 
Einstein radius by 
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M/Mq 



(1) 



where Dl is the distance to the lens, M is the lens mass, 
and we have assumed the source is much further away from 
us than the lens. In practice, the procedure requires the 
Einstein radius to be sufficiently large and the signal-to- 
noise ratio sufficiently high. 

From the presently known catalogues of neutron stars 
and pulsars (see section 4), it appears rare that a galaxy is 
located in the vicinity of a pulsar or a neutron star. For ex- 
ample, in our search results, the smallest angular separation 
between seven ROSAT discovered isolated neutron stars and 
the background galaxies is ~ 50". Although these stars are 
nearby and have relatively high proper motions (hundreds 
of milli-arcsecs per year), such a separation will still require 
a neutron star more than fifty years to overlap with the 
background galaxy even if it moves in the right direction. 

We explore an alternative if we have only one image of 
the background galaxy when the foreground star is already 
well aligned with the background object. If the (unlensed) 
surface brightness of the background galaxy is smooth on the 
scale of the lens distortions (~ 0.1"), it is plausible that the 
observed distortions are due to the lens, from w hich we can 
deriv e the Einstein radius of the lens star (see I Coles et al.l 
2010). In practice, we may have two images where the lens 
moved a detectable distance within a few years, which will 
allow a difference image to be obtained; we return to this 
point in the discussion. 

Let us consider an unlensed image S of the background 
galaxy. It can be expanded in terms of basis functions B" (9) 
as 



S = ^a n B n {6), 



(2) 



where a n, s are the expansion coefficients. Since lensing con- 
serves surface brightness, the lensed image of the galaxy is 



(3) 



where z~t is the position of the lens star and the point <j> in 
the source is given by the lens equation: 



Zt 



(4) 



Zt 



Now we have a model of the lensed image D. Assuming 
the lens position is precisely known, we have one param- 
eter 9e, which we would like to determine, and a set of 
parameters a"'s, which are not of interests. Once we have 
an observed lensed image dij at t when the lens star has a 
position zl, the maximum likelihood is given by 



L(a n ,6 E ) =IJexp[-iff < /(dy - Dij) 2 } 



(5) 



where cry is the noise, ij runs the pixel. We assume the 
measurement errors are Gaussian. Because of this assump- 
tion and the fact that the model D is a linear function of the 
expansion coefficients a n , we can marginalise the coefficients 
a n and only leave one parameter 6 E in the marginalised like- 
lihood (see eq. 8 in lColes et ai]|201Ch . For clarity, we rewrite 
the marginalised likelihood as: 



X 2 = 2\nL{6 E ) = ln|detC|+^P m P n C m?l -^a^4,(6) 



where 



P n — a ij 2 dijL" j 



and 



— 1 \ * —2 j m t n 

^mn / ij ij ij ' 



{-) 



(8) 



The Einstein radius 6 E is involved in the terms L™ through 



(0) 



which are the "lensed" basis functions. Once the noise level 
cry is given from the observed lensed image dij, we can ob- 
tain the best-fit Einstein radius 8 E by minimising the \ 2 ■ 

With this method, even though we do not know the un- 
lensed image, we can still derive the Einstein radius 9 E from 
one observed lensed image. There is, however, an underlying 
assumption that the distortion around the lens is not due to 
any "substructure" in the background galaxy. 



3 SIMULATIONS 

In the section, we explore the feasibility of deriving the Ein- 
stein radii of the pulsars or isolated neutron stars and study 
the corresponding errors. 

First, we produce a simulated lensed image of the back- 
ground g alaxy. We assum e the unlensed surface brightness 
follows a Ide Vaucouleursl (|l948l 'l profile. We would like to 
stress that the specific form of the unlensed image of the 
galaxy is not important as long as it is smooth on a scale a 
few times the Einstein radius of the lens star. In practice, it 
remains to be seen whether this condition is satisfied in real 
images. 

The pair search in Section 4 suggests that the back- 
ground source will be faint and the photons from the sky 
background dominate the noise. Here we consider a spe- 
cific example. We assume the source galaxy has a mag- 
nitude J = 19.86 and the sky surface brightness is J = 
16.75 mag/" 2 , appropriat e for the darkest nights on Mauna 
Kea (|Sanchez et alj|2008f ) . It turns out that in the region we 
consider (0.284" x 0.284") the flux from the sky background 
is 12.5 times that from the source galaxy. So in our simula- 
tion, the noise comes from two parts: one from the source 
itself and one from the sky background. Both follow the 
Poisson distribution. We then obtained the simulated lensed 
image at each pixel dij from this Poisson distribution. The 
mean of the Poisson distribution is the corresponding lensed 
brightness Sy (converted into photon counts) plus photons 
from the sky background. The lensed brightness Sy is given 
by 



S^) = exp(-7.67[V'(^/ J R e ) 2 + R 2 C ] 1/4 ) 



(10) 



with an effective radius R B = 1" and a core radius R c of 
2 pixels (which mimics seeing). The connection between 6ij 
and cj>ij in Eq. (|10[) is given by Eq. Q. The noise cy is taken 

We specialise to the case for the next-generation ex- 
tremely large optical-IR telescopes, in particular the TMT. 
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For the first-light instrument, the InfraRed Imaging Spec- 
trograph (IRIS), the pixel size is 0.004" and the diffraction- 
limited resolution is 0.007" at 1 micron. In our simulation, 
we use the IRIS pixel size as the minimum scale in the choice 
of the scale parameter j3 in the basis functions (see below) . 

We assume the lenses are nearby stars in the Milky 
Way and the background galaxies are far away. For the 
known closest is olated neutron star , the measured distance 
is about 160 pc jPosselt et alj [20091 ). if we assume its mass 
is 1.4M , then 8 E « 0.008". For the pulsar (J1238+21) 
which is crossing in front of a galaxy right now, using a 
mass of 1.4Mq and the estimated distance 1.7 kpc, we ob- 
tain 6 E « 0.0025". Considering TMT has a resolution of 
0.007", we take 8 E = 0.005" and 8 E = 0.01" in our simula- 
tions to produce the lensed images. 

We use shapelets as the basis functions. As we know, 
a faithful shapelets decomposition depends on four parame- 
ters: the scale parameter f3, the maximum order of shapelets 
order -/ V max and the choice of the centroid position of the 
galaxy (|Melchior et alj|2007h . The minimum and maximum 
sizes of fea tures that can b e resolved by the shapelets are re- 
spectively (|RefregieiH2003h min « /3(n tota i + l)~ 1/2 ,# ma x ~ 
/3(itotai + l) 1 ^ 2 , where retotal is the sum of the maximum 
orders for two dimensions (N x ,max and Ny >ma x)- 

In our tests, we find the scale parameter /3 affects the x 
fitting the most. An arbitrary choice of /3 can cause a mis- 
estimation of the coefficients a n . Since the set of a n repre- 
sents the unlensed image, a mis-estimation of a n can lead to 
a mis-representation of the unlensed image. Although a n 's 
are marginalised out, our tests show the different choices of 
P can still affect the behaviour of the effective \ 2 ■ Therefore, 
it is necessary to find an optimal p in advance. In our tests, 
we found the value of optimal f3 changes with the size of the 
pixel and the Einstein radius. In this simulation, we use the 
pixel size of 0.004", and 8e — 0.005". The maximum order 
of the basis function is set to N mBX = 21 for both two di- 
mensions (x and y). T he centre of the image is fixed at the 
centre of the grid box. iRefregieU (|2003h suggested that for a 
given min and max a good choice for /3 is « (0 mia 9 max ) 1/2 . 
Following this, we choose the minimum scale as the pixel size 
(0.004") and the maximum scale as the box size (0.284"), 
and the final chosen scale parameter is /3 = 0.03". 

When the Einstein radius is as small as 0.005", the area 
affected by lensing is small (if we consider 20 times the Ein- 
stein radius, that is 0.1"). If the background galaxy has a 
small gradient in the surface brightness on this scale, the 
lensing effect is too small to be detectable. This is easily 
understood: lensing conserves surface brightness, and a uni- 
form source has no lensing effects. This is the case when the 
lens is far away from the centre of the galaxy. Furthermore, 
the signal-to-noise (S/N) ratio will be low. 

In contrast, near the centre of the galaxy, the surface 
gradient is high, and so even if the Einstein radius of the 
lens star is small, the distortion effect at this place can still 
be strong enough to be detectable. In our simulation, we as- 
sume the lens star is close to the centre of the background 
galaxy. We only use a part of the lensed image which cov- 
ers 71 x 71 pixels. With a pixel size of 0.004", the area is 
0.284" x 0.284". We assume that within this area, the num- 
ber of photons from the galaxy is 3.5 x 10 6 and that from 
the sky background is 12.5 times higher. The centre of the 
galaxy is at (0, 0), and the lens is put at (2.5, —2.5) in units 




Figure 1. The simulated lensed image by a lens star of Einstein 
radius E = 0.005" in the central part of the background galaxy 
which covers an area of 0.284" x 0.284". The effective radius of 
the galaxy is 1". The position of the lens is the centre of the 
circle and the radius of the circle is twice of the Einstein radius. 
The distance between the lens star and the centre of the galaxy 
is 0.014". Assume the photons from sky background noise is 12.5 
times that from the galaxy in this area. The sky background light 
has been subtracted in this image. 



of pixel. This position of the lens corresponds to a distance 
of 0.014". 

The simulated lensed image is shown in Figure [1] The 
distortion within the circle in the image is caused by a fore- 
ground lens star with an Einstein radius 8e = 0.005". The 
centre of the circle is the position of the lens star and the ra- 
dius of the circle is 26 e ■ In this case, the minimum \ 2 value 
is 6873. The degree of freedom is 71 x 71 - (22 x 22 + 1) = 
4556. Then the reduced minimum effective \ 2 is 1-5. Fig- 
ure [2] shows the reduced effective \ 2 curve. The best-fit 
8 E = 0-005^q'qqp2„ . The la error is derived from the scatters 
in 300 Monte Carlo realisations. In Monte Carlo realisations, 
the simulated lensed image dkj are resampled from Poisson 
distributions. So If there are no errors involved in the es- 
timates of distances of the lens star and the source galaxy, 
such an error from 8 E (~ 6%) would translate into an error 
of ~ 12% in the mass estimate of the lens star (8e oc vM). 

Several factors affect the value of the accuracy of the 
best-fit 8e- These are the size of the Einstein radius, the 
total number of photons and the position of the lens. For 
comparison, we have repeated our Monte Carlo simulations 
for several other combinations of parameters. The results 
are shown in Table [1] In this table, the maximum projected 
distance between the lens star and the centre of the galaxy 
is ~ 0.024". It clearly demonstrates that when the Einstein 
radius increases or the lens is closer to the centre of galaxy 
(where the gradient in the surface brightness is high), the 
Einstein radius can be more accurately determined. Further- 
more, when the total number of photons increases, the error 
bar decreases. 
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Figure 2. The reduced effective \ 2 for the case shown in Figure 
[T] The degree of freedom is 4556. The vertical dash line indicates 
the true Einstein radius 8e = 0.005" used in producing the lensed 
image. 



Table 1. ltr errors for different Einstein radii and lens positions. 
The distance in the third row is the angular separation between 
the foreground lens and the centre of the background galaxy. It 
is in units of the Einstein radius. For each simulation, the total 
number of photons from the source galaxy (in the considered re- 
gion) is 3.5 X 10 6 . The photon number from the sky background 
is 12.5 times higher. Be is the true Einstein radius used to pro- 
duce the lensed images. The la errors are given by Monte Carlo 
simulations. 



6 E 


0.005" 


0.005" 


0.01" 


0.01" 


lens 


(5.5,-2.5) 


(2.5,-2.5) 


(5.5,-2.5) 


(2.5,-2.5) 


Distance 


4.8 


2.8 


2.4 


1.4 


ltr error 


+0.0004" 
-0.0003" 


+0.0003" 
-0.0002" 


±0.00020" 


±0.00015" 



4 SEARCHING EXISTING CATALOGUES 

We use the known puls ars in the ATNF Pulsar catalogue^ 
^Manchester et ai"1l2005l ) as our target lens stars. We cross- 
match this catalogue with several galaxy catalogues, in or- 
der to look for existing or promising pulsar-galaxy pairs. 
An existing pair means that the pulsar is already passing 
across the surface of the background galaxy while a promis- 
ing pair means that the pulsar will move across the back- 
ground galaxy within less than one hundred years. The data 
sets include 

(i) ESO archive data from NACO imaging. The search 
radius is 1'. We did not find any pair. 

(ii) 2MASS Extended Source Catalog. We found 10 
pulsar-galaxy pairs with separation of less than 5". Upon 
closer examination, we find that seven of these belong to 
the globular cluster 47 Tuc, and two others belong to the 
globular cluster NGC6544. In these cases, the globular clus- 
ters are the "extended" sources rather than galaxies. Only 



http: / / www.atnf.csiro.au / people / pulsar / psreat / 



in the single remaining pair, the extended source is a faint 
galaxy. The information about the promising pair is listed 
in Table [2] 

(iii) the SDSS catalogue 0. Within a radius of 3", there 
are 4 pairs. The positions of these four pairs are listed in 
Table [3] The best example is shown in Fig. [3] (the second 
row in Table [3j . The other galaxies in the remaining pairs 
all have very low signal-to-noise ratios. Further observations 
are needed to confirm whether they are true galaxies. 

We also searched for possible galaxies surrounding the 
seven thermally emitting radio-quiet iso lated neutron star s 
discovered in the ROSAT all-sky survey (|Motch et al.ll2009l ). 
The galaxy search was limited within the above three 
databases. Among these, J0420. 0-5022 has the closest neigh- 
bour galaxy, ~ 45" away from the st ar. Considering this star 
has a proper motion of 123 mas/yr |Motch et al J 12009(1 . we 
have to wait for more than 100 years for the star to move 
across the galaxy. 



4.1 "lensing" probability and integration time 

The above search suggests that it is rare that a known pulsar 
or a neutron star has a close neighbour galaxy or is pass- 
ing across a background galaxy. Since there will be deeper 
sky surveys, both in the radio (the Square Kilometer Array 
[Sk43), and in the optical (e.g. LSST0, PanSTARRS^), it 
is expected that more and more fainter galaxies and new 
pulsars will be discovered in the future. Promising pairs can 
then be observed with diffraction-limited, extremely large 
optical-IR telescopes. 

In this section, we will estimate, with deep imaging, the 
number of galaxies that can be found in the area of sky swept 
out by the known pulsars within 10 to 15 years. From Table 
[1] we see that when the distance is 4.8#e, we can obtain a 
best-fit 6e from one lensed image with one sigma error bar 
of ~ 8% . In the following, we consider an area affected by 
a lens star is covered by a diameter of 1Q6e- The area swept 
by each pulsar per year is fi x (10 9e), where is the proper 
motion. 

First we estimate the total area for all the pulsars in 
ATNF catalog. In this catalog, there are 241 pulsars with 
known proper motions. Among these pulsars, 70% of them 
has a proper motion within (0, 30) mas/yr. 10 mas/yr is the 
peak value of the proper motion histogram. We then assume 
the remaining pulsars without proper motion measurements 
all have proper motions of 10 mas/yr. With this assumption, 
the total area swept by the pulsars in ATNF catalog within 
15 years by 10 E is 26.51 arcs 2 for 6 E = 0.01" and 13.25 
arcs 2 for 9 E = 0.005". 

We then use the GalaxyCount prograrrjf| 
|Ellis fc Bland-Hawthornl |2007j l to estimate the num- 
ber of galaxies given the magnitude limit. We convert 
the K-band number counts to the J-band since there is 
no J-band option in their program. If TMT can reach 



www.sdss.org 

www . skatelescope . org 

www.lsst.org 

http: / / pan-starrs . ifa. hawaii . cdu/ public / 

The catalog we used was updated on Jul 25, 2012. 

http: / /www. aao.gov.au/astro/GalaxyCount 
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Table 2. List of one pair found in 2MASS extended Source Catalog cross-matched with the catalogue of [Manchest er et al. | (I2005T) . 



5 



Pulsar Name RA Dec Galaxy(NED Obj Name) RA Dec Separation Magnitude 

(hms) (dms) (hms) (dms) (arcsecs) (J band) 

B1900-06 19:03:37.939 -06:32:21.94 2MASX J19033714-0632191 19:03:37.1 -06:32:19 3.59 13.161 



Table 3. List of four pairs found by cross-matching the SDSS catalog and the catalog of iManchester et al. | J2005I) . 



Pulsar Name 



RA 
(hms) 



Dec 
(dms) 



Galaxy(SDSS obsID) 



RA Dec Separation Magnitude 

(hms) (dms) (arcsecs) (r band) 



J0927+23 
J1238+21 
J2317+1439 
B0820+02 



09:27 
12:38 
23:17 
08:23 



37.000 23:46:59.999 587741491435078548 09:27:37.057 23:47:00.675 1.987 22.5 

23.177 21:52:11.075 587742013285794672 12:38:23.203 21:52:11.258 0.414 21.5 

09.237 14:39:31.219 587730774421078668 23:17:09.230 14:39:31.182 0.104 22.8 

09.768 01:59:12.412 588010358525657901 08:23:09.787 01:59:12.581 0.331 23.2 



the magnitude 24.48 in the J-band, and we use a colour 
J — K = 1.0 from UKIRT (|La Barbera et al.ll201oT ). then 
we get a magnitude limit of K = 23.48. The program 
gives 115 galaxies per sq. arcmin with K ^ 23.48. For 
6e = 0.01", within the area, 0.85 galaxies will be found. 
For 9e = 0.005", the rate is reduced by one half to about 
0.42 galaxies. 

It is interesting to estimate the required integration 
time for such a program. The average number of photons is 
about 15 photons s _1 m~ 2 at J = 20 assuming a telescope 
throughput of about 30%. In Table(4] we list the integration 
times required for jtotai — 1.4 x 10 s for five different mag- 
nitudes of the galaxy J=18., 19.496, 20., 20.28, 22. -y total 
is a total number of photons from the galaxy within the 
simulated area 0.284" x 0.284". At faint magnitudes, the 
telescope time would become rather prohibitive. 

We also investigate the observational feasibilities of the 
known pairs listed in the tables. We assume the background 
galaxy has the magnitude 19.50 in the J-banc0 and an effec- 
tive radius of 1", which is similar to the galaxy in the second 
case in Table [3] For this magnitude, the ratio of the light 
from the sky background to the galaxy itself within a central 
area of 0.28x0.28 square arcsecs is 9.0. The distance between 
the centre of the galaxy and the lens is assumed to be 0.105", 
which is close to the third case in Table [3] With a photon 
number of 3.5 x 10 6 from the galaxy in the given central 
area, the lensing effect for an Einstein radius 9e = 0.005" 
is buried in the sky background. But with a higher photon 
number, such as 1.4 x 10 s , which corresponds to a S/N ~ 49 
at the lens position, this method can reproduce the input 
value 6e = 0.005". However, to collect this number of pho- 
tons within the given central area of the galaxy, we need 
about 21 hours for TMT assuming a throughput of 30%. 



5 SUMMARY AND DISCUSSION 

In this work we studied the feasibility of determining the 
mass of isolated neutron stars from their lensing effects on 

9 We use a colour r — J = 1.954 jLa Barbera et al.llioiol) for the 
galaxy. It may have an uncertainty of about 1 mag because we 
do not know the type of this galaxy. 



background galaxies. When a neutron star is passing in front 
of a background galaxy, the surface brightness of the back- 
ground galaxy suffers from distortions. This effect can be 
used to derive the Einstein radius of the lens if we assume 
there are no intrinsic substructures close to the location of 
the lens. 

Our calculations suggest that once the position of the 
lens star is identified precisely, the Einstein radius can be 
derived from only one lensed image of the galaxy given a 
sufficient number of photons. The Einstein radius can be ef- 
fectively derived when the background galaxy has a sharp 
surface brightness gradient near the lens star. We showed 
that the Einstein radius can be derived with 4-8% accuracy 
for a size of 0.005" if we can achieve the required signal-to- 
noise ratio. The resulting mass uncertainty is about 8-16%. 
The main limitation of this method is likely the expensive in- 
tegration time required, even on extremely large telescopes, 
from hours to tens of days. This long integration time may 
be remedied in the future when we have larger pulsar cata- 
logues from SKA, or the overlapping galaxy happens to be 
bright. Such examples may be already seen (see Table 3 and 
Fig. 3). 

We are being somewhat conservative here by taking just 
one image. In fact, for a relative proper motion of 10 mas/yr, 
the pulsar would have moved about 20 mas in two years, five 
IRIS pixels on TMT and 3 times the point spread function. 
One can compare the differences in the two images, and thus 
obtain tighter constraints on the lens mass. We can apply 
our method to the first-epoch image to infer the lens mass, 
and then take a second image to further confirm and verify 
the model. 

If the background source is visible in the radio, then 
an interesting possibility is to observe it with the Square 
Kilometer Array (SKA). In this case, the resolution can be 
of the order of milli-arcsecs, and the dynamical range can 
also be very high. However, one potential problem may be 
the foreground pulsar will be bright in the radio. Similarly, in 
the optical, we assumed the neutron star does not contribute 
any light to the images of the lensed galaxy. However, it 
may not be the case in reality. Therefore, We need to select 
the optimal bands to take images of the galaxy in order to 
minimise the contamination of light from the neutron star. 
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Figure 3. The SDSS image of the galaxy (Obj ID:587742013285794672). Its r-band magnitude is 21.45. At the current time, the pulsar 
J1238+21 is crossing the surface of the galaxy. The pulsar's position in 1999 (J2000): RA 12:38:23.17, Dec +21:52:11.1. There is no 
proper motion information available in the literature and no record shows it is a binary system. The white cross indicates the location 
of the pulsar; the error in the position is about 1 arcsec. 



Table 4. Integration time in hours needed to achieve the signal-to-noise ratio for TMT with 30% throughput for 'Jtotal — 

1.4 x 10 8 . 

Itotal is a total number of photons from the galaxy within the simulated area 0.28" X 0.28". 



magnitude (J band) 18. 19.496 20. 20.28 22. 

Integration time (hours) 5.3 20.8 33.2 42.9 209.02 
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